Transitions between turbulent states in rotating Rayleigh-Benard convection 
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Weakly-rotating turbulent Rayleigh-Benard convection was studied experimentally and numer- 
ically. With increasing rotation and large enough Rayleigh number an abrupt transition from a 
turbulent state with nearly rotation-independent heat transport to another turbulent state with 
enhanced heat transfer is observed at a critical inverse Rossby number l/Ro c ~ 0.4. Whereas for 
1/Ro < l/Roc the strength of the large-scale convection-roll is either enhanced or essentially un- 
modified depending on parameters, its strength is increasingly diminished beyond l/Ro c where it 
competes with Ekman vortices that cause vertical fluid transport and thus heat-transfer enhance- 
ment. 

PACS numbers: 47.27.te,47.32.Ef,47.20.Bp,47.27.ek 



Turbulence evolves either through a sequence of bifur- 
cations, possibly passing through periodic and chaotic 
states [l[ as in Rayleigh-Benard (RB) convection Q when 
the Rayleigh number Ra (to be defined below) is in- 
creased, or through subcritical bifurcations Q as in pipe 
or Couette flow. Once the flow is turbulent, it usually is 
characterized by large random fluctuations in space and 
time and by a loss of temporal and spatial coherence. 
For the turbulent state common wisdom is that the large 
fluctuations assure that the phase-space is always fully 
explored by the dynamics, and that transitions between 
potentially different states that might be explored as a 
control parameter is changed are washed out. 

Contrary to the above, we show that sharp transitions 
between distinct turbulent states can occur in RB con- 
vection 0, Q when the system is rotated about a vertical 
axis at an angular velocity f2. In dimensionless form the 
angular velocity is given by the inverse Rossby number 
1/Ro = 2tt/^(3gA/L. Here L is the height of a cylin- 
drical RB sample, (3 the thermal expansion coefficient, A 
the temperature difference between the bottom and top 
plate, and g the gravitational acceleration. At relatively 
small Ra where the turbulence is not yet fully developed, 
we find that the system evolves smoothly as 1/Ro is in- 
creased. However, when Ra is larger and the turbulent 
state of the non-rotating system is well established, we 
find that sharp transitions between different turbulent 
states occur, with different heat-transfer properties and 
different flow organizations. Similar sharp transitions 
between different states were reported recently for tur- 
bulent flows in liquid sodium [6j], where the increase of 
the magnetic Reynolds number beyond a certain thresh- 
old led to the spontaneous creation of a mean magnetic 
field and where sharp bifurcations between different tur- 
bulent states were observed when a control parameter 



was tuned. One sees that sharp transitions, usually re- 
garded as characteristic of low-dimensional systems, are 
displayed also by the fully developed turbulent flows. 

We present both experimental measurements and di- 
rect numerical simulations (DNS) for a sample with di- 
ameter D equal to L. They cover different but overlap- 
ping parameter ranges and thus complement each other. 
Where they overlap they agree very well. Without or 
with only weak rotation, for this system it is known that 
there are thermal boundary layers (BLs) just below the 
top and above the bottom plate, with a temperature drop 
approximately equal to A /2 across each. The bulk of the 
system contains vigorous fluctuations, and in the time 
average a large-scale circulation (LSC) that consists of a 
single convection roll with up-flow and down-flow oppo- 
site each other and near the side wall. 

The numerical scheme was already described in refs. 
0, H E3]- The experimental apparatus also is well 
documented 0, ELI and we give only a few relevant 
details. The sample cell had D = L = 24.8 cm, with 
plexiglas side walls of thickness 0.32 cm and copper top 
and bottom plates kept at temperatures T t and Tt, re- 
spectively. The fluid was water. The Rayleigh number 
Ra = f3gAL 3 /(Kis) (v and k are the kinematic viscosity 
and the thermal diffusivity, respectively), Prandtl num- 
ber Pr = v I 'k, and Ro were computed from the fluid 
properties at the mean temperature T m = (T t + Tfc)/2. 
The Nusselt number Nu = A e ///A was determined from 
the effective thermal conductivity X e ff = QL/A (Q is 
the heat-current density) and the conductivity A(T m ) 
of the quiescent state. Eight thermistors, labeled k = 
0, . . . , 7, were imbedded in small holes drilled horizon- 
tally from the outside into but not penetrating the side 
wall [Hj]. They were equally spaced around the circum- 
ference at the horizontal mid-plane [z = 0). A second 
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and third set were located at z = —L/4 and z = L/4. 
Since the LSC carried warm (cold) fluid from the bottom 
(top) plate up (down) the side wall, these thermistors de- 
tected the location of the upflow (downflow) of the LSC 
by indicating a relatively high (low) temperature. To de- 
termine the orientation and strength of the LSC, we fit 
the function 

T ftk (z = 0)=T Wi0 + 5 cos(kn/4-e y, fc = 0,...,7 (1) 

separately at each time step, to the eight temperature 
readings Tk (z = 0) obtained from the thermistors at z = 
0. Similarly we obtained 9 t , <5t, and T w _ t for the top level 
at z — L/4. At z = —L/4 only the mean temperature 
T Wt b was used in the current work. 

In Ref. [13] we explored Nu as a function of Ra, Pr, 
and Ro in a large parameter regime, ranging towards 
strong rotation (1/Ro 1) and from small to large Pr. 
Here we focus on Pr « 4 — 7 (typical of water) and weak 
rotation (Ro > 1) to study the transition from the non- 
rotating state at 1 /Ro = towards the rotating case for 
different Ra. 

We start with numerical results for the relatively small 
Ra = 4 x 10 7 which is not accessible with the current ex- 
perimental apparatus because L is too large (and thus 
Ra too small to be accessible with reasonable A). Those 
simulations where done on a grid of 65x 129x 129 nodes in 
the radial, azimuthal and vertical directions, respectively, 
allowing for a sufficient resolution of the small scales both 
inside the bulk of turbulence and in the BLs adjacent to 
the bottom and top plates where the grid-point density 
was enhanced [1, . The small Ra allowed for very long 
runs of 2500 large-eddy turnover-times 2L/U (where U 
is the free-fall velocity U = \//3gAL) and thus excellent 
statistics. Fig.Q]shows the ratio of Nu(fl) in the presence 
of rotation to Nu(il = 0) as function of l/Ro. This ratio 
increases rather smoothly with increasing rotation. This 
increase is thought to be due to the formation of the 
Ekman vortices which align vertically and suck up hot 
(coldjfluid from the lower (upper) BLs (Ekman pump- 
ing) h], 13, 14, 15]. This is supported by the change in 
character of the kinetic BL near the bottom and top walls 
based on the maximum root-mean-square (rms) velocities 
in the azimuthal direction. For 1/Ro < 0.5 the BL thick- 
ness (based on the rms azimuthal velocity) is roughly 
constant or even slightly increases, and for 1/Ro > 0.5 it 
behaves according to Ekman's theory and decreases with 
increasing rotation rate, see the inset in Fig. [1] As the 
turbulence is not yet fully developed, a smooth transi- 
tion between the non-rotating and the rotating turbulent 
states is observed. 

Both numerical and experimental findings are very dif- 
ferent for the larger Ra = 2.73 x 10 8 and Pr = 6.26 where 
the turbulence of the non-rotating system is well devel- 
oped. In Fig.[2]one sees that now there is a critical inverse 
Rossby number l/Ro c ss 0.38 at which the heat-transfer 
enhancement suddenly sets in. For weaker rotation the 
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FIG. 1: The ratio Nu(Q)/Nu(Q = 0) as function of 1/Ro for 
Ra = 4 x 10 7 and Pr = 6.26. Open black squares indicate 
the numerical results. The numerical error is approximately 
0.2% which is indicated by the size of the symbols. Inset: 
The thickness of the kinematic top and bottom BLs based 
on the maximum rms azimuthal (upper symbols: red (black) 
for top (bottom) BL) and radial (lower symbols: blue (green) 
for top (bottom) BL) velocities. The vertical dashed lines in 
both graphs represent l/Ro c and indicates the transition in 
boundary-layer character from Prandtl-Blasius (left) to Ek- 
man (right) behavior. 



data are consistent with no heat-transfer modification 
as compared to the non-rotating case. For this larger 
Ra one sees that experimental and numerical data (now 
based on a resolution of 129 x 257 x 257, see [l(|) agree 
extremely well. In Ref. 13( data from DNS were reported 
on the relative Nusselt number for higher Raylcigh num- 
ber, Ra = 1 x 10 9 and Pr = 6.4, which - as we now 
realize retrospectively - showed a similar transition also 
at l/Ro c w 0.4. 

To characterize the flow field, we numerically calcu- 
lated the rms velocities averaged over horizontal planes 
and over the entire volume, respectively. The maximum 
rms azimuthal and radial velocities near the top and bot- 
tom wall have again been used to define the thickness of 
the kinetic BL, which is shown in the inset of Fig. [2] 
for Ra = 2.73 x 10 8 and Pr = 6.26. The critical in- 
verse Rossby number clearly distinguishes between two 
regimes: one with a constant BL thickness (in agree- 
ment with the presence of the LSC and the Prandtl- 
Blasius BL) and another one with decreasing BL thick- 
ness for 1/Ro > 0.38. The scaling with rotation rate 
is in agreement with Ekman BL theory \ U /L ~ Ro 1 / 2 . 
Further support can be found in Ref. [161 ] where similar 
data for the kinetic BL thickness have been plotted for 
Ra = 1 X 10 9 and Pr = 6.4. For 1/Ro > l/Ro c the nor- 
malized (by the value without rotation) volume-averaged 
vertical velocity fluctuations w rms strongly decrease, in- 
dicating that the LSC becomes weaker, see Fig. [3] The 
decrease in normalized volume averaged vertical velocity 
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FIG. 2: Nu(ft)/Nu(fl = 0) for Ra = 2.73 x 10 8 and Pr = 
6.26. Red solid circles: experimental data (T m = 24°C and 
A = LOOK). Open black squares: numerical results. The 
experimental error coincides approximately with the symbol 
size and the numerical error is approximately 0.5%. Inset: 
Thickness of the kinetic BL. For dashed vertical lines and 
inset: see Fig. [T] 



FIG. 3: The normalized averaged rms vertical velocities w r ms 
for & = 4x 10 7 (left) and Ra = 2.73 x 10 s (right) as func- 
tion of 1/ Ro. The black line indicates the normalized volume 
averaged value of w rms . The red and the blue line indicate 
the normalized horizontally averaged KVms at the edge of the 
thermal BL based on the slope at respectively the lower and 
upper plate. The vertical dashed lines again indicate the po- 
sition of l/Ro c . 



fluctuations coincides with a significant increase of the 
horizontal average at the edge of the thermal BLs indicat- 
ing enhanced Ekman transport (see also insets in Figs. [T] 
and [2|) . These averages provide additional support for 
the mechanism of the sudden transition seen in Nu and 
indicate an abrupt change from a LSC-dominatcd flow 
structure for 1/Ro < l/Ro c to a regime where Ekman 
pumping plays a progressively important role as 1/Ro 
increases. 

Our interpretation for the two regimes is as follows: 
Once the vertical vortices organize so that Ekman pump- 
ing sucks in the detaching plumes from the BLs, those 
plumes are no longer available to feed the LSC which 
consequently diminishes in intensity. A transition be- 
tween the two regimes should occur once the buoyancy 
force, causing the LSC, and the Coriolis force, causing 
Ekman pumping, balance. The ratio of the respective 
velocity scales is the Rossby number. For Ro 3> 1 the 
buoyancy-driven LSC is dominant, whereas for Ro <C 1 
the Coriolis force and thus Ekman pumping is stronger. 
The transition between the two regimes should occur at 
Ro = 0(1), consistent with the observed Ro c ~ 2.6. 

One wonders of course why the transition in between 
the two regimes is sudden (in Nu) for Ra — 2.73 x 10 8 
and less abrupt for the smaller Ra = 4 x 10 7 shown in 
Fig. [T] Our interpretation is that for that relatively low 
Rayleigh number the flow is not yet fully turbulent so 
that spatial coherences [B| exist throughout the cell. For 
1 /Ro < 1 / Ro c the vortical structures have time to form 
and survive for a while in the LSC-induced wind for lower 
Ra, and provide some weak Ekman pumping and thus 
enhanced Nu, while they are easily swept away by the 
much stronger wind of the LSC for higher Ra. In the 



latter case no Nu enhancement is thus observed. 

At higher Ra = 9.0 x 10 9 (where A is larger and tem- 
perature amplitudes can thus more easily be measured) 
and Pr — 4.38 an even more complex situation is re- 
vealed, as seen in Fig. 0] (here DNS is not available be- 
cause it would be too time consuming). We find that 
now Nu(fl) /Nu(Sl = 0) (Fig. fjji), after a slight increase, 
first decreases, but these changes are only a small frac- 
tion of a percent. Then Nu undergoes a sharp transition 
at l/i?o c ,2 = 0.415 (vertical dotted line in Fig. 2£i) and 
beyond it increases due to Ekman pumping. Compari- 
son with Figs. Q] and [2] shows that the transition of Nu is 
not strictly at a constant l/Ro c , but that Ro c depends 
weakly on Ra and/or Pr. 

The LSC amplitudes So and St determined from fits of 
Eq. |T]) to the sidewall-thermometer readings are shown 
in Fig. [4ja as solid symbols. Consistent with the results 
reported in Ref. |l2j, St < So when there is no rota- 
tion (1/Ro = 0). This inequality disappears as 1/Ro 
increases. Both amplitudes first increase by nearly a fac- 
tor of two. At l/Ro c .\ ~ 0.337, where the two amplitudes 
have just become equal to each other, they begin to de- 
crease quite suddenly and remain equal to each other up 
to the largest 1/Ro. The transition at l/i?o Cj i is indi- 
cated by the leftmost vertical dotted line in Figs. and 
c. At that point there also is a transition revealed by the 
vertical temperature difference AT W = 2 x [T Wt b — T Wt t] 
along the side wall as seen in Fig. HJ; which shows AT W /A 
as a function of 1/Ro. Consistent with the initially en- 
hanced LSC amplitudes So and St, these results first show 
a reduction of the thermal gradient as the LSC becomes 
more vigorous, but then reveal an increase due to en- 
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FIG. 4: Results for Ra = 9.0 x 10 9 and Pr = 4.38 (T m = 
40.00°C, A = 16.00 K). (a): Nu(Q,)/Nu(Q = 0) vs. l/_Ro. 
The error bar is smaller than the size of the symbols, (b): 
Solid symbols: time-averaged LSC amplitudes {So}/ A (z = 0, 
circles) and (St) /A (z = L/4, squares) as a function of 1/Ro. 
Open symbols: rms fluctuations about the cosine fit (Eq. [TJ 
to the temperature data, (c): Vertical temperature variation 
AT W /A along the sidewall. (d): Circles: time-averaged nor- 
malized sidewall-temperature profile {[T(9) — T w ]/So) at the 
horizontal midplane for 1/Ro = 1 determined as in JT^]. Solid 
line: cos(0 — 0o). 



hanced plume and/or vortex activity above \/Ro c ,\. 

Also of interest are the rms fluctuations 5T/A = 
([T k {z = 0)-T ftk {z = 0)] 2 ) 1/2 /A about the fit of Eq. Q 
to the temperature measurements at the horizontal mid- 
plane (z = 0), and similarly at z = L/A. They are shown 
as open symbols in Fig. 2)d. These fluctuations begin to 
rise at \/Ro Ci i rather than at \/Ro c .\. Then they soon 
become comparable to 8$ and 5t , suggesting that the LSC 
becomes more and more hidden in a fluctuating environ- 
ment. Nonetheless, remnants of the LSC survive and can 
be found when the fluctuations are averaged away as 
shown in Fig. |4ji. There we see that even for 1/Ro = 1.0 
the time average ((Tk(z = 0) — T Wi q)/Sq) of the deviation 
from the mean temperature T w> q retains a near-perfect 
cosine shape. 

From these measurements we infer that the establish- 
ment of the Ekman-pumping mechanism is a three-stage 



process. First, up to l/Ro c ,i, the time-averaged LSC 
amplitudes, such as (5q) / A, nearly double in value (see 
Fig. 2b) and thereby reduce the vertical thermal gradient 
along the wall (see Fig.0];). Beyond l/Ro Ci \ there is an 
enhanced accumulation of plumes and vortices, which co- 
incides with an increase of the BL thickness near onset as 
shown by the simulations at lower Ra (see insets in Figs.Q] 
and[2]). This accumulation detracts from the driving of 
the LSC but the flow is not yet organized into effective 
Ekman vortices. This organization sets in at \/Ro c .2, 
leads to Ekman pumping, and enhances Nu and reduces 
the strength of the LSC as supported by the volume- 
average of w rms , see Fig. [3] (for lower Ra). This sequence 
of events is altered as Ra (or presumably also Pr) is 
changed, but it is remarkable that for fully developed 
turbulent RB convection sharp bifurcations occur simi- 
lar to those observed in turbulent flows in liquid sodium 
by Q . These resemble the characteristic series of bifur- 
cations well-known from low-dimensional chaos. 
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